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Abstract. We provide a determinantal formula for tau-functions of the KP 
hierarchy in terms of rectangular, constant matrices A, B and C satisfying a 
rank one condition. This result is shown to generalize and unify many previous 
results of different authors on constructions of tau-functions for differential 
and difference integrable systems from square matrices satisfying rank one 
conditions. In particular, it contains as explicit special cases the formula of 
Wilson for tau-functions of rational KP solutions in terms of Calogero-Moser 
Lax matrices as well as our previous formula for KP tau functions in terms of 
almost-intertwining matrices. 



1. Introduction 

In many recent papers by different authors, determinantal formulas have been 
used to transform constant square matrices satisfying a rank one condition into tau- 
functions for integrable systems. In particular, we recall the result of G. Wilson 
that n x n matrices X and Z satisfying the "almost-canonically conjugate" 
condition 

(1) T&vk([X,Z] +7) = 1 

produce tau-functions for rational solutions to the KP hierarchy by the formula 

(2) r(t) = det ^X + ^ iUZ^^j . 

(See also El El HI ■) 

The previous result can be interpreted as a relationship between the KP hierarchy 
and the Calogero-Moser particle system and is therefore similar to the relationship 
between KdV solitons and the Ruijsenaars-Schneider particle system 01 EJ- I n 
that case one finds that the Lax matrices X and Z for this particle system satisfy 
the rank one condition rank (XZ + ZX) — 1 and that the formula 

Coo oo \ 

exp(^ t 2l+1 Z 2i+1 )X exp(^ t 2l+1 Z 2i+1 ) + I 
i=0 i=0 ) 

gives a tau-function for a multi-soliton solution of the KdV hierarchy. 

Both soliton and rational solutions were produced by a formula in our previous 
paper in which it was shown that square matrices X, Y and Z satisfying the 
"almost-intertwining" condition 

(4) rank (XZ - Y X) = 1 
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produce KP tau-functions through the formula 

/ oo 

(5) r(t) = det f X exp( V UZ l ) + exp( 



=1 




Previous results have also demonstrated the usefulness of matrices satisfying 
these rank one equations in constructing solutions to difference equations. In par- 
ticular, it was shown in that for matrices X and Z satisfying QJ, the eigenvalues 
x m f f ne ma t r i x 

(6) X(m) = -r]X ■ (Ai - Z) - mr/(A 2 - Zy 1 • (Ai - Z) 

(where 77 and Ai are arbitrarily selected constants) satisfy the rational nested Bethe 
Ansatz equations 

(7 TT — / 3 — -it- = -1 Vl< j < n. 

t=i ( x 7 ~ x *~ + ^ ( x ? - x k-v) (x? - 4 n+ ) " " 

Moreover, we announced at NEEDS 2001 (an immediate consequence of that 
the formula 

(8) 

r™ = det [X( Cl - Z)\c 2 - Z) m (c 3 - Z) n + (ci - Y)\c 2 - Y) m (c 3 - Y) n W] 
gives a solution to the Hirota Bilinear Difference equation 

/r\\ / \ 771,71 771+1,71+1 / \ 771+1,71 771,71+1 , / \ 771,71+1 771+ 1,71 n 

(9) (c 2 - c 3 )t 1+1 r, - (c x - c 3 )t ( r J+ i + (ci - c 2 )r, t ;+1 = 

when X, Y and Z satisfy J1J (for arbitrary n x n matrix W). 

In addition, two formulas have recently been published for Baker-Akhiezer func- 
tions of the g-KP hierarchy in terms of matrices satisfying g-variants of these "al- 
most" operator identities. In particular, the formula 



det (xzl + xqZ — zX — XZ — I) 
(10) ^' Z) = det(xI-X)det(zI + q Z) " 



XZ 



in the case rank (A" Z — qZX + I) = 1) is found in j^j. Similarly, in JHI one finds 
the formula 



„ s , , , det(qZX - xZ — zX + zxl) 

(11) ip(x,z) = - i-e, 

y ' YK ' ' det(ZX - xZ - zX + zxl) q 



XZ 



for the case rank (XZ — qZX) = 1. In both cases, the interest in these particular 
Baker-Akhiezer functions lies in their bispectrality (cf. [§]). 

Besides the obvious similarities, some explicit connections have been drawn be- 
tween these results. For example, the formula © is easily derived from J3J) (al- 
though this is not the way it was derived in ^U]), formula © is able to construct 
the rational solutions like J2J) as well as soliton solutions, and © can be viewed as 
a special case of . The formulas concerning the g-KP hierarchy should be able to 
be related to the others by the correspondence in |S] (cf. PP) between tau-functions 
and wave functions of the q-KP and KP hierarchies. However, what is clearly lack- 
ing is a general framework in which all of these different formulas appear as special 
cases. It is our goal in this paper to provide such a framework by generalizing the 
formulas to the case of rectangular matrices. 
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2. Main Results 

Let A, C be full rank n x N matrices and let B be an N x N matrix for N > n. 
For convenience, we will assume the non-degeneracy condition det [AC T ] 7^ is 
satisfied. 

Let L = tow A be the subspace in spanned by the rows of A and L 1 - be 
its orthogonal complement with respect to the standard bilinear form (x, y) — 
Sfc=i x iVi- Fi x a basis u\, . . . , ii^-n of L and define U to be the (N — n) x N 
matrix with rows ui, . . . , un~u- As in let g(x) = JZiSa ti xl be a power series 
in x with coefficients that depend on the time variables t = (ti,t 2 , ■ ■ ■) of the KP 
hierarchy. 

Theorem 2.1. If 

(12) rank(AB£/ T ) < 1, 
then 

(13) r™'™ = det [A(al - B) l (c 2 I ~ B) m (c 3 I - B) n C T ] 
is a solution to the Hirota Bilinear Difference Equation jnj and 

(14) r{t) = det 

is a tau-function of the KP hierarchy. 

Proof. It is known (see, e.g. ^21201) that to prove that r(t) is a KP tau-function 
it is sufficient to prove that T/ m ' n := r(t — l^ 1 ] — mfc^ 1 ] — ^[c^ 1 ]) solves the 
HBDE for all values of the parameters. 1 Hence, our method will be to prove the 
second claim above by proving the first. Moreover, it is important to note that 
it is suficient to prove that HBDE is satisfied when all of the discrete "times" 
l,m,n = 0. This is because rj n ' n — Tq'° if f is the tau-function corresponding to 
the same choice of A and B but with a different C (multiplied on the right by the 
transpose of (c±I — B) 1 • •• (03 — B) n ). Since the condition in the claim depends 
only on a property of A and -B, it is therefore sufficient to consider the restricted 
version of the equation 

(02 - c 3 )r(t- [cr 1 ])r(f- N- 1 - [c^ 1 ]) - ( Cl - c 3 )r(t- [c 2 " 1 ])r(f- [q 1 ] - [c, 1 ]) 

+ (c 1 -c 2 )r(t-[c^})T(t-[c^}-[c^}) = 0. 
We will reduce this equation to the following identity, proved in j 1 1 j : 

(15) hi(c 1 )h 2 {c 2 ,c 3 ) - hi(c 2 )h 2 (d,c 3 ) + hi(c 3 )h 2 (ci, c 2 ) = 0, 

where h\{c\) — det[ci - P], h 2 (ci,c 2 ) = det[(ci — P){c 2 — P) + Q] with P,Q n x n 
matrices and rank(Q) < 1. 

First, let V T be any right inverse of A and define G = [V T U T ],B = G~ l BG, 
C T = G- 1 C T and M = M(t) = 

det [M(t)] . Note also that G _1 = 

As usual, the Miwa shift t + c[z] is defined by 

t + c[z] = (ti+cz,t 2 + — ,t 3 + —-■■)■ 



Ae^C T 



[In 

A ' 



e g(v>C 1 . Then AG = [I n 0] and r(t) = 
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The Miwa shift r(t — [c 1 ]) can be computed as follows: 



r(f- 



) = dot [[7 n 0]e s ^e ln ^- c 1 ^C T 



(t) det \c - [I n 0]Be 9{6) C T M(i)- 



and, similarly, 

r(t- [c^ 1 ] - [c^ 1 ]) = Ci n C2 n r{t) det \[I n 0](ci - B)(c 2 - B)t?&6 T M®~ X 

Then it is not hard to check that a left hand side of the bilinear difference Hi- 
rota equation is proportional to the left hand side of ((T5|) . if one defines P = 
[I n 0}Be g ^C T M- 1 and Q = [I n 0}(Bfe 9 ^C J "M" 1 - P 2 . Hence, it suffices to 
show that rank(Q) < 1. Rewrite Q as 

Q = [I n 0]B (i - e 9 ^C T M- 1 [I n 0}) Be 9 ^C J M' 1 



M 

* 



Therefore, 



and notice that e 9(s) C T = 

I -e 9{B) C T M- 1 [I n 0] = 





In-u 



which means that 

[I n 0]B(l- e 9 ^C T M- 1 [I n 0]) = ( [/„ 0}B 



In-t 



[* I 



N-n\ 



But the factor in parentheses is equal to ABU and, by our assumption, has a rank 
less or equal to 1. Therefore, the same is true for Q, which finishes the proof. □ 

Remark: It is worth emphasizing that the rank-one condition i|12|) does not depend 
on a choice of bases in subspaces L and that correspond to the matrix A. 

Since we have constructed a solution to the KP hierarchy, one may be interested 
in knowing where in the Sato-Segal- Wilson grassmannian I17j the correspond- 
ing solutions lie. This can be resolved by studying the associated Bakcr-Akhiezcr 
function, ip(x, z). In particular, one knows that the point W £ Gr is the subspace 
spanned by ip and its ^-derivatives evaluated at x = 

We will show below that there exist polynomials p(z) and q(z) such that p{z)H + C 
W C q^ 1 (z)H+, which shows by definition that W is in the subgrassmannian Gr rat 

0HHj. 

Theorem 2.2. The stationary Baker- Akhiezer function corresponding to the tau- 
function given in <|14l) is 

n( * lf , det(Ae* B (zI-B)C) xz 

(16) ip(x,z) = 



det(Ae xB C) 



As a consequence, we are able to determine that this solution corresponds to a 
point in the subgrassmannian Gr rat of rank one KP solutions with rational spectral 
curves. 

Proof. The formula for ip(t, z) follows from the well known formula for the time- 
dependent wave function |17) : 



il>(t, z) = 



r(t 



it) 



4 



Rank One Conditions for Rectangular Matrices 



Gekhtman/Kasman 



which simplifies when evaluated at t = (a;, 0,0,0, . . .) to the formula above since 
the coefficients of — [x] are precisely the coefficients in the power series expansion 
of log(l — x). 

Another way to identify the subspace W is by its duality with the dual Baker- 
Akhiezer function ip* which by similar arguments to above can be shown to have 
the property that p(z)ip* is nonsingular in z for p{z) = det(z/ — B). Hence, the 
innerproduct of any polynomial in p(z)H + with ip* (computed as the path integral 
of the product around S 1 ) is zero. Consequently, p{z)H + C W and we see that W 
is in Gr rat . □ 

It has frequently been found to be useful in the case of solutions corresponding 
to rational spectral curves [7| |S1 1171 IT%] to identify a solution instead by the finite 
dimensional space of finitely supported distributions in z that annihilate the Baker- 
Akhiczcr function. A consequence of the role of the characteristic polynomial p(z) 
in the proof above is that the finitely supported distributions in z which annihi- 
late z ip(x,z) are supported at the eigenvalues of B with highest derivative taken 
bounded by the algebraic multiplicity of the eigenvalue. 

3. Special Cases 

As a corollary to the main theorem above, we can determine the following gen- 
eralization of our theorem from : 

Theorem 3.1. Let X be an n x (N — n) matrix, Y be an n x n matrix and Z be 
an (N — n) x (N — n) matrix such that 

rank (XZ - Y X) < 1 

then <|14|) is a tau-function of the KP hierarchy where A = [X I n ] , B is the block 
diagonal matrix B = diag[Z, Y] and C is an arbitrary full rank n x (N — n) matrix. 

To see that this is so, note that U in l|12|) can be chosen to be [— I n X T ]. One 
then finds that condition (|12f> reads rank(AZ — YX) < 1. This generalizes (@J, 
since there it was assumed that N = 2n and X is a square matrix. In the latter 
case, one can define C to be C = [I n I n ], which transforms (|14(l into a tau-function 
r(t) = det [Xe^ z ~> + e^ Y ~>] that coincides with ©. 

In fact, combining the discussion of the previous paragraph with the observation 
that the distributions annihilating the Baker-Akhiezer function are supported at 
the eigenvalues of B with order bounded by the algebraic multiplicity confirms the 
conjecture in JI] that the same was true for the eigenvalues of Y and Z. 

Moreover, one can also rederive Wilson's formula (J2J) and the " almost-canonically 
conjugate" rank one condition as a special case of our main result. Consider the 
case when N — 2n, A, C and U arc defined as in above but B is chosen in the form 

Z ° Then ABU T = - ([X, Z] + I n ) and thus, JEJ) coincides with JTJ. 



B 



In Z 

Moreover, in this case e 9 ^ = 



and the tau-function 114|) 



e 9(Z) 

g'{Z)e3<- z ) 

becomes r(t) = det[e 9(z) ] det[X + g'(Z)], which is gauge equivalent to the one in 
©• 

In conclusion, although many of the details are yet to be fully explored, the 
formula we have proven above for the first time allows us to consider many different 
results relating rank one conditions and tau- functions in a unified context. We 
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plan to address questions of the relationship between the geometry of the space 
of matrices we utilize and the geometry of the grassmannian, reductions to finite 
dimensional Hamiltonian (particle) systems, and the case in which A, B and C are 
taken to be infinite dimensional operators in a future paper. 
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